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ABSTRACT 

The existence of planets born in environments highly perturbed by a stellar companion represents 
a major challenge to the paradigm of planet formation. In numerical simulations, the presence of a 
close binary companion stirs up the relative velocity between planetesimals, which is fundamental in 
determining the balance between accretion and erosion. However, the recent discovery of circumbinary 
planets by Kepler establishes that planet formation in binary systems is clearly viable. We perform 
N-body simulations of planetesimals embedded in a protoplanetary disk, where planetesimal phasing 
is frustrated by the presence of stochastic torques, modeling the expected perturbations of turbulence 
driven by the magnetorotational instability (MRI). We examine perturbation amplitudes relevant to 
dead zones in the midplane (conducive to planet formation in single stars) , and find that planetesimal 
accretion can be inhibited even in the outer disk (4-10 AU) far from the central binary, a location 
previously thought to be a plausible starting point for the formation of circumbinary planets. 

Subject headings: Planets and satellites: formation, Planets and satellites: dynamical evolution and 
stability 



1. INTRODUCTION 

Among the numerous discoveries of the Kepler 
mission, the detection of circumbinary (CB) plan- 
ets Kepler 16 through the multiple sy s tem Kepler 
47 dPovle et all 1201 it [Welsh et all I2012t lOrosz et al.l 
l2012all bl: iSchwamb et al.l I2012T ) has propelled a renewed 
theoretical effort in explaining how such planets would 
be assembled in a binary environment. In presence 
of a close binary companion, protoplanetary disks can 
become a rather hostile planetary nursery in both cir- 
cumstellar and circumbinary configurations. The dif- 
ficulties encountered by the standard core accretion 
paradigm in the binary environment are several, includ- 
ing truncation, mass loss an d relatively fas t dispersal of 
disks in close binaries (e.g. lDuchenell2010t iKraus et al.l 
12012ft . possible vaporization of grains in dynamically 
excited disks ([Nelson! I2000f) a nd impaired planetesi- 
mal growth into protoplanets ([Marzari fe Scholll 1 2000 1: 
Moriwaki fe Nakagawall2004t iThebault et alJl200l boOft 
Scholl et all l2007t IPaardekooper et al.l 120081 : iThebaultl 
2011). 

The latter "planetesimal bottleneck" is a robust conse- 
quence of the interplay between the gravitational pertur- 
bations of the stellar companion (which stirs the plan- 
etesimal disk and acts to raise eccentricities) and the 
aerodynamic drag from a putative protoplanetary disk. 
Since aerodynamic drag tends to both damp planetesi- 
mal eccentricities and align planetesimal orbits in a size- 
dependent fashion, planetesimals of different sizes will 
collide on different phases, resulting in high collisional 
speeds which lead to destruc t ive (r ather than accreting) 
events. Recently, iMeschiaril (|2012fl (hereafter M12) in- 
vestigated planetesimal accretion in the Kepler-16 sys- 
tem, using A-body simulations (coupled with drag from 
a static background disk) which track planetesimal col- 
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lisions throughout a range of semi-major axes over 10 5 
years. The census of planetesimal collisions indicated 
that regions inside 4 AU (« 20a#, where as is the semi- 
major axis of the central binary) were dominated by de- 
structive events, and therefore hostile to planet forma- 
tion. Consequently, we posited that Kepler 16 could have 
plausibly assembled outside the forbidden region and 
subsequently migrated inwards through tidal interaction 
with the protoplanetary disk, later stopping close to the 
inner edge of the disk. This scenario is supported by 
the evolu tion of planetary cores i n hyd rodynamical sim- 
ulations (jPierens fe Nelson! [2007ll2008D . which also sug- 
gested Jupiter-mass planets in CB configurations should 
be rare (in accordan c e with the observed Kepler sample). 
IPaardekooper et al.1 (|2012f ) investigated planet formation 
in the Kepler 16, 34 and 35 systems using a similar ap- 
proach, additionally including self-consistent planetesi- 
mal formation and destruction. They reached analogous 
conclusions, and asserted that in- situ formation was un- 
likely, even under the most favorable conditions. 

The approach of Ml 2, similarly to previous investi- 
gations, neglected several physical ingredients for the 
sake of computational expediency. Indeed, fully self- 
consistent simulations which include the hydrodynamical 
response of the protoplanetary disk have shown that ad- 
ditional oscillations in eccentricity and longitude of peri- 
center of the planetesimals might be introduced from the 
development of bulk eccentricity and spiral perturbations 
in the disk. However, their magnitude might depend 
somewha t on the details of the hydrodynamical simula- 
tion (e.g. IPaardekooper et al.|[2008t iMarzari et al.ll2Cjr2t 
IMuller fe KlevH2012D . 

In this Letter, we consider magnetohydrodynami- 
cal turbulence driven by the magnetorotational in- 
stability (MRI) as an additional source of perturba- 
tions on the planetesimal dis k. MRI-driven turbu- 
lence (|Balbus fc Hawlevl I1991D is thought to be the 
likely source of anomalous viscosity in protoplanetary 
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disks (e.g. lArmitagel 1 1 9981). and influence ho w plan- 
etesimal formation (e.g. Johansen et al.l |2007[). plan- 



etesimal accretion (e.g . lOgihara et al.l 12007; Id a et all 
120081 : INelson fe Gressell2010h. and planetary mig ration 
(e.g. iNelson fe Papaloizoul 120031 lLaughlin et alJ I2004t 
INelson fe Papaloizoul 12004 iBaruteau fe Linl 120101 ) pro- 
ceed. Previous studies simulating the dynamics of plan- 
etesimals embedded in turbulent disks in single-star en- 
vironments showed that in fully MRI-active disks, the 
velocity dispersion of the planetesimals is significantly 
raised by gravitational perturbations induced by density 
fluctuations. For typical disk parameters and nominal 
turbulence strength, km-sized planetesimals might be in 
a highly erosive regime. Such vigorous turbulence, how- 
ever, might not be appropriate to the midplane of real- 
istic protoplanetary disks, which is thought to be domi- 
nated b y a "dead zone " with near-laminar flow (|Gammid 
11996ft . iGressel et all (|201lD (hereafter Gil) presented 
the results of stratified, 3D MHD simulations (includ- 
ing a substantial dead zone) ; planetesimals embedded in 
the disk midplane at 5 AU experienced a significantly 
reduced excitation of their eccentricities (by a factor w 
10-20). Therefore, they concluded dead zone represent 
"safe havens" for the growth of km-sized planetesimals. 

Although the reduced amplitude of the random veloc- 
ities excited by turbulent fluctuations is potentially con- 
ducive to planetesimal accretion in the single-star envi- 
ronment, the situation is more complicated in the binary 
environment. In the latter case, a very precise alignment 
between the planetesimal orbits is crucial to attaining 
low encounter speeds despite the substantial eccentricity 
of the planetesimal orbits. Random kicks diffuse plan- 
etesimals out of alignment; although aerodynamic drag 
would attempt to restore alignment, it will do so on a 
timescale that is size-dependent, once again differentially 
phasing planetesimals of different sizes. Finally, high- 
frequency radial oscillations of the planetesimal eccen- 
tricity (as would be caused by stochastic torques) would 
potentially lead to further orbital crossing, leading to 
high encounter velocities. Therefore, we anticipate that 
planet formation might be strongly perturbed, or even 
inhibited, despite the smaller turbulent amplitudes ap- 
propriate to the midplane dead zone. We couple the 
./V-body code of Ml 2 with a numerical model that ap- 
proximates the stochastic torquing arising from MRI tur- 
bulence. Our analysis shows that this additional source 
of perturbations is potentially damaging to planet for- 
mation in the outer disk, which was previously thought 
to be relatively protected from destructive impacts and 
therefore a plausible loc ation for core assembly (Ml 2; 
iPaardekooper et aI1l2012l ). 

The plan of this Letter is as follows. In $2j we briefly 
discuss our numerical model for the gravitational torques 
arising from MHD turbulence. In [J3]we show the results 
of our simulations, and discuss them in the context of 
planet formation in 2) 

2. NUMERICAL SETUP 

In this Letter, we study planet formation by sampling 
planetesimal collision events between 4 and 10 AU, and 
compare them to velocity thresholds corresponding to 
destructive impacts. We consider the orbital elements of 
the Kepler-16 system as our prototypical configuration, 
and to facilitate comparison with previous investigations 
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Fig. 1. — Sample growth plot of the eccentricity dispersion for 
models A, Bl-4 (grey lines, from top to bottom). The dashed 
lines indicate the best-fit eccentricity dispersion measured in the 
full MHD simulations of Gil, models A, Dl and D2 (from top to 
bottom). 



(M12; Paardekoopcr ct al. 2012). We refer the reader to 
M12 for a description of the numerical code, initial setup 
and velocity thresholds. 

2.1. Turbulent model 



We use the analytical prescription of lLaughlin et al.l 
(2004) (hereafter L04) to model torques arising from 
density fluctuations in the disk; we include corrections 
to the formulation from lOgihara et all (120071) (also used 
in llda et a l. 2008; Barutca u fc Linll2010h . and introduce 
small modifications that reduce its computational cost. 
In this model, density fluctuations are forced by a poten- 
tial $, a sum of turbulent m-fold modes 

$ = jr 2 il 2 2J £,iR(r, r c ,i, ai) T(U) cos(mcp-(p c> i-Cl c>i ti) 

i 

(!) 

where 7 sets the overall turbulent amplitude, r c .i, (p c j 
and tt c .i are the radial center, phase and angular velocity 
of the mode (picked randomly in the disk), m is sampled 
from a lograndom distribution between 1 and 6, <Ji = 
Tc,i / 4m is the radial extent of the mode, £j is sampled 
from a Gaussian distribution of unit variance, U is the 
lifetime of the mode normalized by a timescale Atj and 
R and T are two Gaussian-like functions centered around 
r Ci i and 0.5 Ati, respe ctively. Each mode ha s a limited 
lifetime Atj; following lBaruteau fc Linl (|2010l) . we reduce 
the lifetimes by a factor of 10 from the prescription of L04 
in order to better match the autocorrelation timescale of 
3D MHD simulations. 

The actual force on the planetesimals will arise from 
the gravitational force of the density fluctuations induced 
by Equation ([1]). In order to proceed without the full 
hydrodynamical machinery, L04 used a WKB analysis 
to derive the following scaling for the RMS torque on a 
planetesimal of mass M. p v. 



t t = C 7 r 4 f2 2 £ 



Mp 



(2) 



where M.* is the mass of the central object (in our case, 
the total binary mass) and S g is the unperturbed sur- 
face density of the protoplanetary disk. While L04 also 
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derives an approximate value for the constant C, the 
heuristic nature of the derivation (whic h is equivalent 
to a dimensional analysis, as derived by Uohnson et al.l 
2006) suggests a better approach would be to fit the am- 
plitude of the perturbations 7 to the results of MHD 
simulations. 

To compare with the results of Gil, we integrated the 
orbits of an ensemble of 100 planetesimals, started at 
5 AUs with zero eccentricity and random initial phase, 
subject to the stochastic torques described above and no 
aerodynamic drag and computed the diffusion of their 
orbital elements. We consider several models, differ- 
ing only by the turbulent amplitude 7: model A (7 = 
2.5 x 10~ 3 ) represents a fiducial turbulent amplitude ap- 
propriate to fully MRI-active disks, while models Bl-4 
(7 = 2.5 x 10~ 4 , 10~ 4 , 5 x 10~ 5 , 2.5 x 10~ 5 , respectively) 
have reduced turbulent amplitudes appropriate to the 
midplane of dead zones. Since a scales as -f 2 (where a 
is the usual vi scosity parameter in th e Shakura-Sunyaev 
prescription; IBaruteau fc Lin! [20101 ). the values chosen 
span two orders of magnitude of turbulent viscosity. Fi- 
nally, for comparison with M12, we also ran simulations 
with no turbulence (model C). 

Figure [T] shows the growth of the dispersion of the plan- 
etesimal eccentricity er(e) as a function of time; as ex- 
pected from a random-walk process, <r(e) oc i 1 / 2 . Models 
Bl and B2 bracket well the growth of the eccentricity 
dispersion seen by Gil for models with dead zones; we 
include the reduced values of models B3 and B4 to test 
the robustness of our results. 

We remark that it is likely that the interaction be- 
tween the planetesimals and the background disk is more 
complicated than the model presented in this Letter; for 
instance, we ex pect the disk t o be endowed with some ec- 
centricity (e.g. (Marzari et al.ll200l iMarzari et alJl2012f ). 
However, for the sake of simplicity and to highlight the 
role of turbulence as an additional factor in dephasing 
planetesimal orbits, we decided to take the orthogonal 
approach of ignoring the self-consistent evolution of the 
disk. In this picture, we assume a circular disk as a max- 
imally accretion- friendly starting point. 

3. SIMULATIONS 
3.1. Diffusion of eccentricity and dephasing 

To visualize the effect of the turbulent fluctuation on 
a population of planetesimal, we first integrated the tra- 
jectories of a swarm of 10,000 planetesimals uniformly 
distributed between 4 and 8 AUs for four different lev- 
els of turbulence: no turbulence, model A (active MRI- 
disk), model Bl (nominal dead zone turbulence fitting 
the results of Gil) and model B4 (turbulence amplitude 
reduced by a factor of 10). 

Figure [2] shows the eccentricity and longitude of peri- 
center (w) of the planetesimals. Compared to the run 
without turbulence, both models with levels of turbu- 
lence appropriate to a dead zone add a considerable 
amount of noise on top of the e and n7-profiles secularly 
imposed by the central binary. Since the orbits of neigh- 
boring planetesimals will not be as collimated, impact 
velocities will be raised everywhere in the disk. There- 
fore, we expect that the fraction of accreting impacts will 
be decreased even when very low levels of turbulence are 
present. 
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Fig. 2. — Eccentricity and longitude of pericenter w — mg as 
a function of semi-major axis for 10,000 planetesimals distributed 
between 4 and 8 AU, for runs with no turbulence (black dots), 
model A (grey dots, active-MRI disk), model Bl (red dots, dead 
zone) and model B4 (yellow dots, dead zone turbulence amplitude 
reduced by a factor of 10). 

For the model representing a fully MM active disk, ec- 
centricities are raised to very high values (e«5xl0~ 2 ) 
and the longitude of pericenter is completely random- 
ized; therefore, the collision velocity can be directly es- 
timated as Av « evkep > 400 m/s, comfortably above 
any velocity threshold resu lting in destructive impacts 
([Stewart fc Lei nhardt 2009) . Indeed, in this case the tur- 
bulent torques completely overwhelm the secular forcing 
of the central binary, resulting in high impact speeds that 
are consistent with the results of simulations of planetesi- 
mal d ynamics in fully turbulent disks around single stars 
(e.g. llda et all2008l:INelson fc Gressell[2010l) . 

3.2. Collision statistics 

We subsequently ran full simulations with collision de- 
tection between 4 and 10 AU, for each of our models. 
Given the computational overhead of calculating and up- 
dating the turbulent forcing at each time step, we chose 
to instead integrate the trajectories of a smaller number 
of planetesimals (1,000) concentrated in annuli centered 
around 4, 6, 8 and 10 AUs. This approach has the ad- 
vantage of reducing the running time of our simulations, 
while simultaneously improving our collision statistics by 
increasing the impact rate. 

Figure [3] shows the fraction of accreting impacts as 
a function of semi-major axis for each of our models. 
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Fig. 3. — Fraction of accreting impacts as a function of semi- 
major axis (from lighter to darker: models B1-B4 and C.) Each 
bar represents the percentage of unperturbed impacts; the fraction 
of perturbed impacts is additionally represented by a solid line. 

For the sake of completeness, we further classify impacts 
as "unperturbed" (Av < v esc , where v esc is the escape 
velocity) or "perturbed" accretion (v esc < Av < v ero ); 
for the former, the impact velocity is low enough to allow 
runaway growth to proceed, therefore allowing the rapid 
formation of oligarchs. We generously deem a radial bin 
as a favorable location for planet formation when the 
fraction of accreting impacts (unperturbed or perturbed) 
is larger than 50%. 

For models Bl and B2 (with turbulent amplitudes 
which best fit the results of the 3D MHD simulations of 
Gil), we find that the percentage of accreting impacts 
decreases dramatically. This is consistent with the noisi- 
ness of the pericenter phasing and the high-frequency jit- 
ter around the damped value. These short-term oscilla- 
tions are not damped efficiently by the aerodynamic drag 
at these large distances, since the drag torque decreases 
steeply with distance (tq <x a~ 3 ~^, for planetesimals at 
the equilibrium forced eccentricity) while the turbulent 
torque in our model decreases more slowly (tt oc a 1 ^' 3 ), 
where we take £ oc r"' 3 (/3 = 1.75 for the minimum-mass 
solar nebula model considered in this paper). 

For models B3 and B4 (with reduced turbulent am- 
plitudes) we find that, despite the fact that substantial 
oscillations in eccentricity and longitude of pericenter are 
still induced, the resulting impact velocities have crossed 
the accreting threshold. The median impact speed for 
model B4 is ~ 6 m/s, which is in the erosive regime for 
1-km planetesimals but allow accretion for 10-km plan- 
etesimals. 

Finally, in accordance with the simulations of M12, 
we find that planet formation can proceed undisturbed 
outside 4 AU when turbulent forcing is switched off. 

4. DISCUSSION 

In this Letter we have investigated planetesimal accre- 
tion in the outer parts of a circumbinary disk (4-10 AU). 
Our simulations indicate that the stochastic forcing of 
turbulent perturbations will frustrate planetesimal phas- 
ing and raise eccentricities, inhibiting planet formation 
even quite far from the central binary. This result is ro- 
bust for levels of turbulence that match those observed 
in realistic MHD simulations. In runs modeling disks 
fully invaded by MRI turbulence, planetesimal phasing 



is completely destroyed and planetesimals will collide at 
speeds that are much higher than fiducial erosive ve- 
locities. We then ran simulations fitting the turbulent 
amplitude corresponding to the stratified disks endowed 
with a large dead zone obtained in Gil, resulting in a 
reduction by a factor « 10 — 20 in the turbulent torque. 
While the reduced stochastic perturbations would be suf- 
ficiently small to allow accretion in a single-star environ- 
ment, in our circumbinary configuration the fraction of 
impacts resulting in accreting events was greatly reduced 
(< 1 - 30% between 4 and 10 AU for models Bl and 
B2, respectively) due to the increased eccentricity jitter 
and dephasing. Therefore, it is possible to inhibit planet 
formation by two different mechanisms: by high eccen- 
tricities, differential phasing, and perturbations caused 
by the self-consistent reaction of the gas disk at small 
radii, and by stochastic turbulent torques at large radii. 
This is problematic, since if that picture is accurate, it 
is not feasible to form the observed planetary census at 
any realistic distance from the binary. The circumbinary 
environment is less robust than single stars with regards 
to dynamical perturbations induced by the disk (regard- 
less of their origins) ; a further reduced level of turbulence 
may be required to form planets at all. Indeed, we found 
that further reduction of the turbulent amplitude by a 
factor of 2-5 was necessary to make planet formation vi- 
able again in the outer disk. 

Several uncertainties in our model may still provide 
some room for allowing planet formation despite the per- 
turbing effects of turbulence. Our analytic prescription, 
while attempting to match MHD simulations, only pro- 
vides a "0-th" order description of gravitational pertur- 
bations induced by torques. For instance, we remark that 
we tuned the amplitude of turbulent torques based on 
the evolution of swarm of planetesimals in a single MHD 
simulation, and rescaled the turbulent amplitude at each 
radial location according to the numerical prescription 
of Section [5] however, the output of these models will 
depend on the assumed magnetic field, the ionization 
level and possibly the resolution of the shearing box. 
The detailed radial dependence and even extent of the 
dead zone r emain uncertain and depend on the assumed 
model (e.g. lMatsumura fc Pudritzl 120051: iTerqueml 120081: 
iFlaig et al\\20ll ; if the dead zone extends to only a few 
AUs, then fully active turbulence may still play a role 
in the outer disk and even shut off planetesimal accre- 
tion completely). Effectively, using the torque scaling of 
Equation [2] implicitly assumes that the turbulent torque 
is well described by the fully turbulent model of L04, 
attenuated by a constant factor at each radial location. 
Only self-consistent (including the time-dependent po- 
tential of the binary), computationally expensive global 
simulations could, in principle, inform the model pre- 
sented in this Letter. However, since we have considered 
various values of 7 for each independent annulus, we are 
reasonably covering a number of possibilities for the ra- 
dial dependence of the turbulence amplitude. Our results 
at a given radius should only be considered suggestive of 
further challenges to planet formation. 

Finally, we remark that, similarly to Ml 2, we have cho- 
sen to study planetesimal growth through mutual colli- 
sions in the 1-10 km planetesimal regime. A more sophis- 
ticated approach following the evolution of the size dis- 
tribution of the planetesimals (adding a significant com- 
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putational overhead) may be warranted; in the latter 
approach, coevolution with a disk of small dust might 
make accretion possible closer to t he binary (assuming 
a hig h dust accretion efficiency; iPaardekooper et al.1 
2012). Even in this setup, some amount of migration 
is needed to bring the fully formed core (or embryos) 
to the current location. The assumption of a static 
gas background is not realistic, especially in the inner 
disk; in that case, self-consistent gas dynamics will act 
to further increase impac t velocities in the inner disk 
(jPaard ckoopcr et al.l 120081 ). We plan to incorporate the 
full hydrodynamical evolution of the gas disk in a follow- 
up paper using the SPH module of our code. 

Our simulations indicate that planetesimal accretion 
in the 1-10 km range will be inhibited everywhere in 
the disk, and could only proceed if the strength of the 
turbulent torques is reduced from our fiducial value, or 
the initial planetesimal population is comprised of big- 
ger objects. The latter scenario is particularly appeal- 
ing, as it could allow planet formation to proceed in 
other highly dynamically disturbed environments (e.g. 



lThebaultj[20H . Indeed, our fiducial model for dead zone 
turbulence becomes accretion-friendly for a planetesimal 
size spectrum spanning 10-100 km size. Recent simula- 
tions of planetesimal formation in weakly turbulent disks 
show that massive bound clumps are formed rapidly from 
meter-sized boulders withi n pressure bumps, for typica l 
nebula parameters (e.g. Uohansen et al.l [20071 1201 IT) . 
Such bound clumps will likely result in planetesimals 
with sizes comparable to at least a substantial fraction 
of the dwarf planet Ceres. We suggest that a primor- 
dial population of large planetesimals might be crucial to 
proceed with planet formation in highly perturbed envi- 
ronments. However, further simulations will be required 
to assess whether the formation of large clumps is robust 
to the dynamical perturbations of a binary companion. 



S.M. acknowledges insightful discussions with G. 
Laughlin and N. Haghighipour, a useful critique from 
the anonymous referee, as well as support from the W. 
J. McDonald Postdoctoral Fellowship. 



REFERENCES 



Armitage, P. J. 1998, ApJ, 501, L189 
Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 
Baruteau, C, & Lin, D. N. C. 2010, ApJ, 709, 759 
Cuzzi, J. N., Hogan, R. C, & Shariff, K. 2008, ApJ, 687, 1432 
Doyle, L. R., Carter, J. A., Fabrycky, D. C, Slawson, R. W., 
Howell, S. B., Winn, J. N., Orosz, J. A., Prsa, A., Welsh, 
W. F., Quinn, S. N., Latham, D., Torres, C, Buchhave, L. A., 
Marcy, G. W., Fortney, J. J., Shporer, A., Ford, E. B., 
Lissauer, J. J., Ragozzine, D., Rucker, M., Batalha, N., 
Jenkins, J. M., Borucki, W. J., Koch, D., Middour, C. K., Hall, 
J. R., McCauliff, S., Fanelli, M. N., Quintana, E. V., Holman, 
M. J., Caldwell, D. A., Still, M., Stefanik, R. P., Brown, W. R., 
Esqucrdo, G. A., Tang, S., Furesz, G., Geary, J. C, Berlind, P., 
Calkins, M. L., Short, D. R., Steffen, J. H., Sasselov, D., 
Dunham, E. W., Cochran, W. D., Boss, A., Haas, M. R., 
Buzasi, D., & Fischer, D. 2011, Science, 333, 1602 
Duchgne, G. 2010, ApJ, 709, L114 

Flaig, M., Ruoff, P., Kley, W., & Kissmann, R. 2012, MNRAS, 

420, 2419 
Gammie, C. F. 1996, ApJ, 457, 355 

Gressel, O., Nelson, R. P., & Turner, N. J. 2011, MNRAS, 415, 
3291 

Ida, S., Guillot, T., & Morbidelli, A. 2008, ApJ, 686, 1292 
Johansen, A., Klahr, H., & Henning, T. 2011, A&A, 529, A62 
Johansen, A., Oishi, J. S., Mac Low, M.-M., Klahr, H., Henning, 

T., & Youdin, A. 2007, Nature, 448, 1022 
Johnson, E. T., Goodman, J., & Menou, K. 2006, ApJ, 647, 1413 
Kraus, A. L., Ireland, M. J., Hillenbrand, L. A., & Martinache, F. 

2012, ApJ, 745, 19 
Laughlin, G., Steinacker, A., & Adams, F. C. 2004, ApJ, 608, 489 
Matsumura, S., & Pudritz, R. E. 2005, ApJ, 618, L137 
Marzari, F., Thebault, P., & Scholl, H. 2008, ApJ, 681, 1599 
Marzari, F., Baruteau, C, Scholl, H., & Thebault, P. 2012, A&A, 

539, A98 

Marzari, F., & Scholl, H. 2000, ApJ, 543, 328 
Meschiari, S. 2012, ApJ, 752, 71 

Moriwaki, K., & Nakagawa, Y. 2004, ApJ, 609, 1065 

Miiller, T. W. A., & Kley, W. 2012, A&A, 539, A18 

Nelson, A. F. 2000, ApJ, 537, L65 

Nelson, R. P., & Gressel, O. 2010, MNRAS, 409, 639 

Nelson, R. P., & Papaloizou, J. C. B. 2003, MNRAS, 339, 993 

— . 2004, MNRAS, 350, 849 

Ogihara, M., Ida, S., & Morbidelli, A. 2007, Icarus, 188, 522 
Orosz, J. A., Welsh, W. F., Carter, J. A., Brugamyer, E., 
Buchhave, L. A., Cochran, W. D., Endl, M., Ford, E. B., 
MacQueen, P., Short, D. R., Torres, G, Windmiller, G, Agol, 
E., Barclay, T., Caldwell, D. A., Clarke, B. D., Doyle, L. R., 
Fabrycky, D. C, Geary, J. C, Haghighipour, N., Holman, 
M. J., Ibrahim, K. A., Jenkins, J. M., Kinemuchi, K., Li, J., 
Lissauer, J. J., Prsa, A., Ragozzine, D., Shporer, A., Still, M., 
& Wade, R. A. 2012a, ArXiv e-prints 



Orosz, J. A., Welsh, W. F., Carter, J. A., Fabrycky, D. C, 
Cochran, W. D., Endl, M., Ford, E. B., Haghighipour, N., 
MacQueen, P. J., Mazeh, T., Sanchis-Ojeda, R., Short, D. R., 
Torres, G., Agol, E., Buchhave, L. A., Doyle, L. R., Isaacson, 
H., Lissauer, J. J., Marcy, G. W., Shporer, A., Windmiller, G., 
Barclay, T., Boss, A. P., Clarke, B. D., Fortney, J., Geary, 
J. C, Holman, M. J., Huber, D., Jenkins, J. M., Kinemuchi, 
K., Kruse, E., Ragozzine, D., Sasselov, D., Still, M., 
Tenenbaum, P., Uddin, K., Winn, J. N., Koch, D. G, & 
Borucki, W. J. 2012b, ArXiv e-prints 

Paardekooper, S.-J., Leinhardt, Z. M., Thebault, P., & Baruteau, 

C. 2012, ArXiv e-prints 

Paardekooper, S.-J., Thebault, P., & Mellema, G. 2008, Monthly 

Notices of the Royal Astronomical Society, 386, 973 
Pierens, A., & Nelson, R. P. 2007, A&A, 472, 993 
— . 2008, A&A, 483, 633 

Scholl, H., Marzari, F., & Thebault, P. 2007, MNRAS, 380, 1119 
Schwa mb, M. E., O rosz, J. A., Carter, J. A., et al. 2012, 

larXiv: 1210.36121 
Stewart, S. T., & Leinhardt, Z. M. 2009, ApJ, 691, L133 
Terquem, C. E. J. M. L. J. 2008, ApJ, 689, 532 
Thebault, P. 2011, Celestial Mechanics and Dynamical 

Astronomy, 111, 29 
Thebault, P., Marzari, F., & Scholl, H. 2006, Icarus, 183, 193 
Thebault, P., Marzari, F., Scholl, H., Turrini, D., & Barbieri, M. 

2004, A&A, 427, 1097 
Welsh, W. F., Orosz, J. A., Carter, J. A., Fabrycky, D. C, Ford, 
E. B., Lissauer, J. J., Prsa, A., Quinn, S. N., Ragozzine, D., 
Short, D. R., Torres, G, Winn, J. N., Doyle, L. R., Barclay, T., 
Batalha, N., Bloemen, S., Brugamyer, E., Buchhave, L. A., 
Caldwell, C, Caldwell, D. A., Christiansen, J. L., Ciardi, 

D. R., Cochran, W. D., Endl, M., Fortney, J. J., Gautier, III, 
T. N., Gilliland, R. L., Haas, M. R., Hall, J. R., Holman, M. J., 
Howard, A. W., Howell, S. B., Isaacson, H., Jenkins, J. M., 
Klaus, T. C, Latham, D. W., Li, J., Marcy, G. W., Mazeh, T., 
Quintana, E. V., Robertson, P., Shporer, A., Steffen, J. H., 
Windmiller, G., Koch, D. G, & Borucki, W. J. 2012, Nature, 
481, 475 



